Scaling of the magnetic response in doped antiferromagnets 
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A theory of the anomalous ui/T scaling of the dynamic magnetic response in cuprates at low doping is 
presented. It is based on the memory function representation of the dynamical spin suceptibility in a doped 
antiferromagnet where the damping of the collective mode is constant and large, whereas the equal-time spin 
correlations saturate at low T. Exact diagonalization results within the t-J model are shown to support assump- 
tions. Consequences, both for the scaling function and the normalization amplitude, are well in agreement with 
neutron scattering results. 
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Magnetic properties of cuprates, as they evolve by dop- 
ing the reference antiferromagnetic (AFM) insulator and lead 
to a high-T c superconductor, have been so far a subject of 
intensive experimental and theoretical investigations. One 
of the puzzles awaiting proper theroretical explanation is 
the scaling behavior of the magnetic response observed in 
cuprates, mostly in the regime of low doping jlj]. It has been 
first found by the inelastic neutron scattering experiments in 
La 2 _ K Sr x Cu0 4 (LSCO) at low doping x = 0.04 H that 
q-integrated spin susceptibility Xi( w ) follows a simple uni- 
versal behavior in terms of the scaling variable ui/T. Simi- 
lar scaling has been observed also in La2- as Ba :r Cu04 J^t, in 
YBaCu 3 6+x (YBCO) with x = 0.5, 0.6 Jl, and even more 
pronounced in Zn-substituted YBCO fl. 

In particular, experiments on cuprates at low doping indi- 
cate that one can represent results for the local susceptibility 
in a broad range of u> and T as x'Li 10 ^) = I(ui)f(oj/T) 
where I(ui) — x'[ (w, T = 0). The scaling function should ap- 
proach f(x—* oo ) = 1 and the simplest form invoked in the 
analysis is f(x) ~ (2/7r)tan _1 (Ar). It seems, however, that 
the function is not universal for all cases, i.e., A ~ 1.0 — 2.2 
varies between YBCO results \^ and LSCO, whereby in 
the latter case corrections to the simplest form reveal an even 
better agreement |2]. At the same time it is found that the in- 
verse AFM correlation length k = l/£, as extracted from the 
q-dependent x"(q, w )> saturates at low u> and T. The largest 
response is at the AFM wavevector Q = (tt, tt) and as a con- 
sequence of the scaling the peak in x"(Q, u>) should move 
downward with decreasing T, this being in fact established in 
YBCO 1 6 ] . It should be noted that also NMR relaxation exper- 
iments test and confirm the ui/T scaling of at 10 ~* 
10. 

Such aw/T scaling is inconsistent with the concept of usual 
Fermi liquid. This has been recognized quite early and the 
concept of the ' marginal' Fermi liquid has been introduced 1 8 ] 
to explain scaling of the magnetic response as well as of other 
anomalous electronic properties. One appealing explanation 
still considered is the vicinity of the quantum critical point 
0. However, the latter should in general require also a criti- 
cal variation of k(u>, T), as indeed observed in LSCO near the 
optimum doping [10]. A random-phase-approximation treat- 



ment of x(q, Hill yields a relaxation rate T oc l/£ 2 , and 
the scaling form could be reproduced provided that the cor- 
relation length is critical, i.e., £ oc T -1 ' 2 , which is not con- 
sistent with experiments |; lj|2|]. On the other hand, numerical 
investigations of the twodimensional t-J model confirm the 
scaling of Xz,( w ) 11211 . although results are restricted to rather 
high T as compared to experiments. 

In the following we will argue that the anomalous u>/T scal- 
ing of the magnetic response can be understood as a conse- 
quence of few simple ingredients which appear to be valid 
for doped AFM in the normal state: a) the collective mode is 
strongly overdamped, whereby the damping is nearly to- and 
T- independent at low u>, and b) there is no long-range spin 
order at low T, so that static spin correlations saturate with a 
finite £. It will be shown that these prerequisites are sufficient 
to reproduce several experimental findings for xl{w). 

Within the memory function approach 1 13] the dynamical 
spin susceptibility Xq( w ) = can be expressed in 

the form 



X ^ UJ) u 2 + wM q (w) - ujI 



(1) 



suitable for the analysis of the magnetic response, as manifest 
in underdoped AFM [14]. u q represents the frequency of a 
collective mode provided that the mode damping is small, i.e., 
7 q ~ Mq(wq) < u>q. For 7 q > w q the mode is overdamped. 
The advantage of the form Q is that it can fullfil basic sum 
rules even for an approximate M q '. Thermodynamic quanti- 
ties entering Eq. Q can be expressed as 



(2) 



where x q = Xq( w — 0) is the static susceptibility. 

?7 q is closely related to the spin stiffness and can be ex- 
pressed in terms of the static correlation functions, and is ex- 
pected to be weakly q-dependant for q ~ Q. Static x q (or w q ) 
remains to be determined, even for known Af q (w). Instead of 
directly evaluating x% being quite a sensitive quantity, we 
rather fix it by the sum rule 



duj Cth — X'^LO) = (S^) = C^, 



(3) 
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given in terms of equal time correlations, which are expected 
to be less T-dependent. C q are bound by a local constraint 
(l/iV)^ q C q = (1 — c/t)/4, where c/j is an effective hole 
doping. 

Let us define now our central assumptions. We first take 
that static correlations follow the standard Lorentzian form, 



i.e., C q = C/(k 2 + q 2 ) |2], where q = q Q, although 



our results are not very sensitive to the explicit form of C q (at 
fixed k). k is taken as a T-independent constant, at least on 
approaching low T. Such an assumption for n is consistent 
with the neutron scattering data for weakly doped LSCO [2] 
and YBCO flEll, as well as with results for the t-J model 
at finite doping 1 15]. It furthermore indicates that the system 
remains paramagnetic with finite AFM £ down to the lowest 
T, as well as the absence of any ordered ground state. 

Less plausible is the second assumption that the damping is 
also constant, M q (cj) ~ 7, i.e., (roughly) independent of oj, q 
and T, or at least not critically dependent on these variables. 
We can give several arguments in favor of this simple choice. 
Recently the present authors 1 14] studied the spin dynamics 
within the t-J model, 
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(4) 



with the nearest neighbor £y = t and next-nearest-neighbor 
hopping tij = t 1 . It has been shown that the dominant contri- 
bution to the damping M^(u>) in a doped system comes from 
the decay of spin fluctuations into fermionic electron-hole ex- 
citations 0- If the fermionic excitations in a doped system 
behave as in a Fermi liquid, and the Fermi surface crosses the 
AFM zone boundary, the damping in the normal state is es- 
sentially constant at low uj, and also weakly dependent on T 
and q ~ Q. This is clearly very different from an undoped 
AFM where one expects vanishing M q (tj q ) for q — > Q and 

r^oli. 

In order to support the simplification of constant 7, we 
present in the following numerical results for the t-J model, 
obtained via the finite-T Lanczos method (FTLM) 1 12J for 
a system of N = 20 sites on a square lattice with peri- 
odic boundary conditions. The model is analysed for the 
parameter J/t = 0.3 as appropriate for cuprates (note also 
the relevant value t ~ 400 meV), and within the regime 
of low hole doping, Ch = Nh/N < 0.15. Note that re- 
sults within the FTLM have macroscopic relevance for high 
enough T, while at T < Tj s they become influenced by 
finite-size effects. As a criterion for Tf s we use the thermo- 
dynamic sum Z(T) = Tr exp(-(iJ - E Q )/T) and the re- 
quirement Z(Tf s ) = Z* ^> 1 111 211 . In the cases discussed 
here Tf s ~ 0.1 1 at intermediate doping, i.e., Tf s ~ 400 K in 
terms of cuprate parameters fl2ll . Within the FTLM we calcu- 
late directly x q ( w )- Since ?7 q and w q are given as frequency 
moments of x'q(u), it is then easy to extract also the damping 
function A/ q (o;) via Eq.([0. 

In Fig. 1 we present results both for x q ( w ) an d f° r Mq 
for fixed doping c h = 2/20 and T = 0.15 t > T fs . In the 



analysis a smoothing e = 0.07 t is used for convenience. One 
can conclude that in the presented case we are clearly deal- 
ing with overdamped spin dynamics for all presented q. In 
spite of widely different x q ( w ) the damping function M q '(cj) 
is nearly constant in a broad range of uj < t and almost inde- 
pendent of q. For this particular we estimate k = 0.7 so 
that the span of q goes beyond q > k. 
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Figure 1: (a) Spin susceptibility Xqi^) within the t-J model for 
doping Ch — 2/20 and T — 0.15 t, for different q, b) damping 
function M^(ui) for the same parameters. Inset shows nonequivalent 
q for the lattice of N = 20. 

Fig. 1 confirms that it is meaningful to extract 7 q = 
Mq(ui — 0), which we present in Fig. 2 for q = Q and 
various doping Ch = Nh/N < 0.15 as a function of T. As 
expected, our results in an undoped system are consistent with 
vanishing 7q(T — > 0) whereas for c/, > they lead to a fi- 
nite extrapolated value 7q(T — > 0). Note that the slopes of 
Jjq(T) I dT in Fig. 2 are quite similar for all presented c^. 
This can be interpreted as a signature that the damping is a 
sum of the spin-exchange contribution and the fermionic con- 
tribution 1 14]. Only the spin-exchange term is active in an un- 
doped system and apparently it adds to the fermionic damp- 
ing, the latter dominating the T — > behavior in a doped 
system. It should be also noted that the characteristic (satura- 
tion) scale for the dominant spin-exchange damping is T ~ J 
which is far above the T investigated in experiments. Hence, 
for the T window of interest our results confirm that in a doped 
system the simplification of constant 7 is sensible. On the 
other hand, it is evident from Fig. 2 that 7 increases with dop- 
ing, becoming very large 7 ~ t on approaching the 'optimum' 
doping Ch ~ 0.15. 
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Figure 2: Low-frequency damping 7q as a function of T for various 
doping c h = Nh/N. 



Let us consider the consequences of proposed simplifica- 
tions taking both 7 and 77 as constants. The dynamical suscep- 
tibility now takes the resonance form 



(u 2 



(5) 



which we have to investigate together with the sum rule, 
Eq.(|3j. At given q there are several regimes with respect 
to the values of T, 7 and w q . Most evident is the situa- 
tion for 7, T <C ujq with an underdamped mode with a fre- 
quency (xJ q ~ r]/(2Cq) = a(q 2 + k 2 ) where a = r//2C. 
In such a case both x'q(o;) as well as the local susceptibil- 
ity x'[ ( w ) = (V-^O X) q Xq( w ) show more or less (depending 
on 7) pronounced gap for uj < uiq ~ an 2 . This regime 
clearly does not exhibit the desired oj/T scaling. On the other 
hand, even in a weakly doped AFM at low T with k«1 one 
should enter such an underdamped regime for the collective 
mode with q ^> k. Still, in this case for w q > 7 one would 
expect that the dispersion becomes that of AFM paramagnons 
with w q oc q. This indicates that a Lorentzian form for C q 
presumably is not appropriate for such a regime and should 
be modified for g>K. 

Experiments on cuprates as well as numerical results for the 
t-J model (as apparent in Figs. 1,2), however, show that in the 
normal state the collective mode is always overdamped in the 
vicinity of q = Q, i.e., cjq < 7. Now, one gets a simple 
Lorentzian for low u> < uj a , 



7(^ 2 + r q )' 



r = _a 



(6) 



Let us first discuss low T — > 0. In this case the l.h.s. of 
Eq.(|3j can be explicitly integrated, and for w q < 7 we get 
C q ~ (277/777) ln(7/w q ). The relevant quantity is the peak 
frequency uj p = Tq(T — * 0). We see that the crucial param- 
eter is 



C = C 7 r 7 /(2r ? K 2 ), uj v 



7e 



-K 



(7) 



which exponentially renormalizes u> p . Since C is fixed by the 
total sum rule, i.e., C ~ (1 — Cft)7r/(21n(7r/«;)) ~ 0(1) and 
77 ~ 0.6 1 11411 at low doping, £ is effectively governed by the 
ratio 7/K 2 . Our results for the t-J model, as presented above 
as well as the analysis of experiments on cuprates, indicate 
that generally £ ^> 1. 

A nontrivial quantity which is the consequence of the pre- 
sented T = analysis is the local 0) = /(w) directly 
related to the measured 'normalization' function In 
order to evaluate the latter we first find for each q the ap- 
propriate ujq satisfying the sum rule and then integrate 
over q. Results for I(u>) at various ( are presented in Fig. 3. 
For convenience we fix 7 = 0.2 t, which appears to corre- 
spond (see Fig. 2) to low doping c/, ~ 0.05, close to dop- 
ing in cuprates with observed scaling behavior. We note that 
the range ( = 2 — 8 presented in Fig. 3 corresponds to 
K = 0.35 — 0.19. We see from Fig. 3 that the behavior for 
all ( is qualitatively similar at high ui while the difference is 
mainly in the position of ui p where I(uj) is maximum. 
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Figure 3: Local susceptibility x[(u),T = 0) for different f. The 
inset shows the comparison with (intensity scaled) experimental data 
for 'normalization' function I(uS) in LSCO 1 2] and Zn-doped YBCO 
1 5]. The vertical scales are adjusted since the experimental scales are 
not in absolute units. 



and T q < w q . An overdamped form as in Eq.© has been 
frequently invoked in the analysis of the magnetic response 
I2I I1 ill in the normal state of cuprates. However, without the 
knowledge of w q Eq.0 is not sufficient to analyse the relation 
with the sum rule, Eq.(|3j- 



We can make a direct comparison with experiments on 
cuprates at low doping which reveal a nontrivial I(u>), as pre- 
sented in the inset of Fig. 3. We first note that data for LSCO 
at x — 0.04 1 2] and Zn-substitued YBCO |5] are quite similar. 
Both indicate a steep increase of I(u>) below uj = 10 meV and 
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no sign of saturation even at lj = 2 meV. In terms of our anal- 
ysis this means that ( > 1. The comparison of I(uj) with our 
results at £ = 8 (any £ ^> 1 would be in fact quite satisfac- 
tory) reveals very good agreement . The difference seems to 
appear at larger uj > 20 meV which could be again due to our 
Lorentzian form of C q . Namely, taking for g > k C q a 1/q 
would lead to a flat ~ const., as in an ordered AFM 

where only transverse fluctuations - magnons - contribute. 

We next discuss the behavior at T > 0. It is evident 
that for T > uip the temperature dependence of Tq (T) (or 
u>q(T)) becomes crucial. In order to satisfy the sum rule (0 
it follows Tq(T) oc T which is the origin of the u>/T scal- 
ing. In Fig. 4 we present the 'scaling' function J(uj/T) = 
x'/(w,T)/xi(w,0) for various T and chosen ( = 8. Re- 
sults confirm that indeed f{u/T) is nearly universal in a very 



broad range of T, i.e., between T 



and T 



7- 



We 



show in Fig. 4 for comparison also experimental scaling func- 
tion for Zn-substituted YBCO |5] which generally fits our 
results very well. It is also evident that at least at lower 
T < 0.05 t our scaling function can be closely represented 
by f(x) = (2/7r)atan(Aa;) with A ~ 1.2. 

There is still some dependence of f(x) on parameter £. A 
general tendency is that at larger ( > 1 we observe the sat- 
uration at somewhat smaller uj/T ~ 1, i.e., appropriate A 
increases. On the other hand, for decreasing £ — > 1 we get 
f(x — > 0) > and the saturation moves to somewhat higher 
x. 
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oretical confirmation. However, both requirements are inti- 
mately related since they are consistent with a paramagnetic 
liquid, with fermionic excitations dominating 1ow-cj, low-T 
behavior. This picture is picture is supported by ARPES ex- 
periments revealing well pronounced quasiparticle excitations 
even in a weakly doped LSCO I17II . 

In the presented picture the broad validity of the scaling is 
due to large £ ^> 1, i.e., the AFM q = Q collective modes 
are heavily overdamped even at low T. This is consistent 
both with neutron scattering results in cuprates as well as with 
available numerical results within the t-J model. Note, how- 
ever, that our results should remain equally valid for the case 
of strong-coupling, i.e., U ^> t effective single band Hub- 
bard model, which in the low energy sector reduces to the t-J 
model with J = At 2 /U . Our scenario for the uj/T scaling 
differs from a quantum-critical one [9] in spite of a similar 
behavior of Tq oc T, since k does not scale in the same way. 
It should be pointed out, however, that k as deduced, e.g., 
from Xq or from Xq( w ) at fixed uo is significantly reduced, 
i.e., k < k at low T. 

There are still some open questions. The theory predicts the 
existence of the crossover temperature T ~ lo v below which 
the scaling would cease to exist, and the response would ap- 
proach T = 0). Such a saturation has so far not been 
reported for weakly doped LSCO and Zn-substitued YBCO, 
which do not exhibit other phase instabilities at low T. One 
should, however, not forget a possible influence of disorder, 
since the same region of the phase diagram is often associated 
with the spin glass character Since our assumptions ap- 
pear to remain valid also at higher (up to optimum) doping in 
the normal state, we can speculate on a possible validity of 
the same scenario in this regime as well, provided that other 
instabilities are absent (e.g., superconductivity, stripe order- 
ing). The indication for the latter are the NMR T\ relaxation 
results 1 7] showing the same scaling in LSCO for T > T c up 
to x = 0.15. 

Authors acknowledge the support of the Ministry of Edu- 
cation, Science and Sport of Slovenia. 
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Figure 4: Scaling function f(u)/T) for £ = 8 for different T (in 
units of t). For comparison data for Zn-substitued YBCO are also 
plotted, as measured for different energies w, taken from Ref.| 5]. 

In conclusion, we presented a theory giving an explanation 
for the anomalous uj/T scaling behavior in magnetic response 
of doped AFM. It is based on two key assumptions: a) the sat- 
uration of static spin correlations and of the correlation length 
£ at low T and, b) on the constant damping 7 of the collec- 
tive mode. While the former is supported by experimental 
data [|2l|4j|5|], the latter follows from numerical analysis on 
finite clusters within the t-J model and deserves further the- 
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